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Helioseismological sound-speed profiles severely constrain possible deviations from standard solar 
models, allowing us to derive new limits on anomalous solar energy losses by the Primakoff emis- 
sion of axions. For an axion-photon coupling g al ^ 5 x 10 -10 GeV -1 the solar model is almost 
indistinguishable from the standard case, while g al > 10 x 10~ 10 GeV" 1 is probably excluded, cor- 
responding to an axion luminosity of about 0.20 L© . This constraint on g a -y is much weaker than the 
well-known globular-cluster limit, but about a factor of 3 more restrictive than previous solar limits. 
Our result is primarily of interest to the large number of current or proposed search experiments 
for solar axions because our limit defines the maximum g al for which it is self-consistent to use a 
standard solar model to calculate the axion luminosity. 

PACS numbers: 14.80. Mz, 96.60.Ly 



I. INTRODUCTION 

Stars are powerful sources for low-mass weakly inter- 
acting particles. The backreaction of this energy- loss 
channel on the properties and evolution of stars has been 
extensively used to constrain nonstandard neutrino cou- 
plings and the properties of axions and other hypothet- 
ical particles |Q. While the Sun no doubt is the star 
best known to us, it does not provide the most restric- 
tive limits on new particle properties. Still, it remains 
of interest what the Sun can tell us about nonstandard 
particle-physics assumptions. 

Much has improved in our knowledge of the Sun over 
the past decade because its interior properties have be- 
come accessible by helioseismological methods and by a 
host of solar neutrino experiments. It used to be that 
a solar model was fixed by adjusting its presolar helium 
abundance and the mixing-length parameter such as to 
reproduce the observed luminosity and radius at an age 
of about 4.6 x 10 9 years. In this way even rather dra- 
matic modifications could be accommodated. For exam- 
ple, a nonstandard energy-loss channel in the form of 
axion emission could have been essentially as large as 
the solar photon luminosity so that the present-day Sun 
could have been close to its main-sequence turn-off 0. 
Naturally, one expects that the precise helioseismological 
sound-speed profiles which have become available over 
the past few years will provide far more severe restric- 
tions on possible present-day solar models. 

As a specific example for a nonstandard energy-loss 
channel we consider the Primakoff conversion of axions in 
the Coulomb fields of charged particles, 7 + Ze — > Ze + a 
(Fig. [j]). This case is of particular interest because one 
can search for the solar axion flux by the reverse pro- 
cess where an axion oscillates into an X-ray in a long 
dipole magnet which has been oriented toward the Sun, 



the so-called "helioscope" method ||[|. An alternative 
approach uses the Primakoff backconversion in a germa- 
nium crystal where one can achieve a large enhancement 
in some directions in analogy to Bragg diffraction ||. 
In either case it is of interest to know the maximum 
axion-photon coupling strength which is compatible with 
well-established solar properties. For example, an early 
experimental search for solar axions j(| was not sensi- 
tive enough to detect these particles even if the coupling- 
strength had been so large as to push the present-day 
Sun close to its main-sequence turn-off. Put another 
way, the limit found in this early experiment was not self- 
consistent in that it required a larger axion emission from 
the Sun than is compatible with its age. A surprisingly 
large number of current (7||] and proposed ||[l0| solar 
axion search experiments has recently emerged. This ex- 
perimental activity motivates us to examine solar axion 
limits in the light of helioseismological information. 

In Sec. II we describe our standard solar model and 
compare it with recent helioseismological data. In Sec. Ill 
we construct models of the present-day Sun including 
axionic energy losses and derive a new solar limit on the 
axion-photon coupling constant. In Sec. IV we discuss 
our results in the light of other astrophysical limits and 
in the context of solar axion search experiments. 




FIG. 1. Primakoff production of axions in the Sun. 
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II. SOLAR AND SEISMIC MODEL 

The solar model used in this paper was computed with 
the GArching SOlar Model (GARSOM) code, which has 
been described in its numerical and basic physical details 
elsewhere [0. It was shown |l2j that it agrees in its 
gross features with other contemporary up-to-date stan- 
dard solar models. In particular, it uses the latest OPAL 
opacities and equation of state jl^,[l4j and takes into ac- 
count particle diffusion of hydrogen, helium and a num- 
ber of heavier elements (e.g. C, N, O). We usually follow 
the evolution beginning from the pre-main sequence to an 
age of 4.57Gyr ("best model"). To derive the differential 
effect of axion emission on solar models it is, however, suf- 
ficient to calculate models including only H/He-diffusion 
starting from the zero-age main sequence (ZAMS). 

The main difference to other comparable solar models 
is our treatment of the atmosphere and the outermost su- 
peradiabatic convective layers down to an optical depth 
of 1000, where the stratification of a two-dimensional hy- 
drodynamical simulation Jl5| is used. In the deeper layers 
convection is treated according to Ref. |l6| even though 
the temperature gradients are almost perfectly adiabatic 
and thus independent of the convection theory. 

Since the first publication of our solar model Jll| the 
code has been revised. Apart from some minor details 
the following numerical treatments have been changed: 
particle diffusion and the nuclear network are now solved 
simultaneously in the same system of equations instead of 
following a sequence of burning- mixing-burning between 
two models of successive age. The opacity interpolation 
(a two-dimensional bi-rational spline with one free pa- 
rameter) now resembles closely a standard cubic-spline 
interpolation, because we found that the strong spline 
damping we previously used to prevent unphysical oscil- 
lations of the interpolation function degraded the agree- 
ment of our models with seismic ones |17|. Furthermore, 
the reaction rate for 3 He( 3 He, 2p) 4 Hc JL8| and the solar 
radius (R Q = 6.95508 x 10 10 cm fj||) were updated. 

Theoretical solar models are compared with seismic 
ones which are constructed by the so-called inversion 
method from measured p-mode frequencies. In Fig. || we 
show the difference in sound speed of our best solar model 
(solid line) compared to the seismic model from |2(]] . The 
dashed line shows a solar model which was calculated 
from the ZAMS including only H/He-diffusion (gio = 0). 
The latter is used in the following as reference model 
to illustrate the effect of axion emission. The shaded 
area shows the uncertainty of the derived profile of the 
sound speed (c s ) according to Ref. pi]] . They were quot- 
ing the errors in the quadratic isothermal sound speed. 
To get the respective errors in c s we neglect the uncer- 
tainties from the adiabatic index Ti, as these contribute 
barely to the total error ||^] . There are three sources for 
the total uncertainty: (i) the errors of the measured fre- 
quencies; (ii) the dependence of the final seismic model 
on the starting model; and finally, (iii) the uncertainties 



in the regularization procedure of the inversion method. 
The three uncertainties were determined in a conserva- 
tive way: each error interval was doubled and then their 
absolute values were added ^J. This was done because 
the parameter study may not have been exhaustive. The 
measured frequencies being very precise, the main uncer- 
tainties arise from the inversion method itself. Previous 
studies by ]22| ] found an uncertainty of approximately 
5 x 10~ 4 for 0.2 < r/R Q < 0.8 for each of the three 
error sources quoted above. Applying the same proce- 
dure for determining the conservative errors gives a value 
of 3 x 10 -3 slightly bigger than the uncertainties found 
by pl]| . p3] | remark that seismic models do not really 
constrain solar model for t/Rq < 0.05. 

The deviation of our standard solar model from the 
seismic model by |2C]] (solid line in Fig. ||) remains al- 
most everywhere within the error range. The deviation 
is very similar to that of other comparable solar models 
in the literature. In particular, we also find a large devia- 
tion immediately below the convective zone which might 
be indicative for either a slight error in the opacities or 
for overly effective helium diffusion out of the convective 
envelope (24]. 
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FIG. 2. Difference in sound-speed profiles of present-day 
solar models compared to the seismic model. The shaded are 
reflects the uncertainties in the infered sound speed of the 
seismic model. The solid line shows our best solar model, 
the dashed line the reference model (gio = 0) used for our 
analysis. 

In the first row of Table [0] we list the expected count- 
ing rates for the gallium and chlorine experiments as well 
as the flux of 8 B neutrinos at the (Super)Kamiokande de- 
tector. Assuming flavor conversion of two neutrino types 
as the solution of the solar neutrino puzzle we show the 
parameter space in Fig. |5| for which the measured neu- 
trino fluxes of all experiments are reproduced. The up- 
per panel shows the MSW-solutions ^5|, the lower one 
the solution for vacuum oscillations (shaded areas). In 
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the calculation of the MSW-solutions we include the av- 
erage earth- regeneration effect Q , and for the vacuum- 
oscillations the annual Sun-Earth distance variation. Our 
95% confidence level regions agree with those in the lit- 
erature. 

We have verified by test calculations that including 
metal diffusion and the new value for 5i7(0) of the 
7 Be(p,7) 8 B reaction of 0.019 keV barn @ instead of 
0.0224 keV barn j|| leads to almost the same parame- 
ter regions in the Am 2 -sin 2 20-plane as for the reference 
solar model. 



III. SOLAR MODELS WITH AXION LOSSES 

We now calculate self-consistent solar models where we 
include energy losses by axion emission from zero age to 
the present-day Sun including only H/He-diffusion. Ax- 
ions [ p9| generically have a two-photon coupling vertex in 
full analogy to neutral pions. The interaction Lagrangian 
can be written in the form 

£a7 = 9a~i& ■ Ea, (1) 

where g aj is a constant with the dimension (energy) -1 , 
E and B are the electric and magnetic fields, and a is 
the axion field. In keeping with previous works we will 
always use the dimensionless parameter 

g 10 = . 9a7 /10- 10 GeV- 1 (2) 

to characterize the coupling strength. The energy-loss 
rate as a function of gio is discussed in Appendix A. 

In Table | we summarize the characteristics of the so- 
lar models which include axion losses for several values 
of gio- We show the present-day axion luminosity L ai 



TABLE I. Solar models with axion losses. 
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TABLE II. Neutrino detection rates. 
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the presolar helium abundance Y, and the helium abun- 
dance, density, and temperature at the solar center. In 
Fig. |3| we show the deviation of the sound speed profiles 
of the axionic solar models from a model without axion 
loss (gio = 0). This deviation shows a local maximum 
at around r 0.1 i?©, a region where both the mod- 
els are most sensitive to axion emission and the seismic 
models are well-constrained. The fractional deviation 
I Acg/cs.olmax at this local maximum is listed in Table |[ 
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FIG. 3. Difference in sound-speed profiles of present-day 
solar models with axion losses compared to the reference 
model in the sense (Reference— Model)/Reference. Different 
line types correspond to different values of the axion-photon 
coupling constant: <?io=4.5 (solid line), 10 (short-dashed), 15 
(dash-dotted), 20 (dash-dot-dot-dotted). The shaded area is 
the same as in Fig. ^. 

In Fig. U we show the helium profiles of our axionic 
solar models. The depth of the convective envelope is 
affected very little by the axion loss — all models are 
within the value predicted from helioseismology of 0.710- 
0.716 Rq j3C| ]. On the other hand, the present surface he- 
lium content depends on the initial one, reduced by diffu- 
sion in the course of the evolution by about 0.03. Axionic 
models with gio > 10 have a surface helium abundance 
significantly lower than 0.238, the smallest value allowed 
from helioseismology [^l| . 

The additional energy sink also affects the present- 
day core structure (Table |) and thus the solar neutrino 
flux. The expected neutrino rates are given in Table || 
for all current experiments. Consequently, the neutrino- 
oscillation parameters Am 2 and sin 2 29, for which the 
measurements are reproduced, depend on giQ. For sev- 
eral values of 510 we show in Fig. || the 95% confidence 
regions for the MSW and vacuum solutions of the so- 
lar neutrino problem. The two MSW-solutions (usually 
called small and large mixing angle solutions) move to- 
ward each other with increasing axion loss. In the case 
of i?io=15 there remains only one solution. 
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FIG. 4. Helium profiles of present-day solar models with 
different axion losses. The insert shows a magnification of 
the region around the bottom of the convective envelope. The 
line styles correspond to those of Fig. |^. 



To illustrate that for g w > 15 the allowed neutrino 
parameters fit only very poorly the experimental val- 
ues we show in Fig. ^ the predicted event rates or neu- 
trino fluxes for GALLEX/SAGE, Homestake and (Su- 
per )Kamiokande. For each value of gio the expected 
rates for the best-fit small-mixing-angle (SMA), large- 
mixing-angle (LMA) and vacuum (VAC) solution are 
shown. The error bars reflect the theoretical uncertain- 
ties in the predicted rates, the shaded band the mea- 
surement errors pl[ . In particular, the result from (Su- 
per)Kamiokande disfavours any neutrino oscillation so- 
lution for gio > 15. Additionally, the absence of an 
observed day-night effect in Supcr-Kamiokande excludes 
the region between the SMA and LMA solution of a stan- 
dard solar model (<?io=0). Including this result would 
lead to SMA and LMA solutions for gio > 10 which re- 
produce the measured rates even worse. 



IV. DISCUSSION AND SUMMARY 

We have calculated a series of self-consistent solar mod- 
els with varying amounts of axionic energy losses. For an 
axion-photon coupling gio < 4.5, corresponding to an ax- 
ion luminosity L a < 0.04 Lq, the small modification of 
the solar model is well within the uncertainties of all cur- 
rent solar observables. On the other hand, for 510 > 10, 
corresponding to L a > 0.20 Lq, the modifications are so 
significant that it is unlikely that they can be compen- 
sated by uncertainties of standard-model input param- 
eters such as the opacities, nuclear fusion rates, age or 
metal abundance. 

Our best solar model deviates from the seismic model 
at most by ss 0.15%, just within the generous unccrtain- 
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FIG. 5. Neutrino mixing parameters in a two-flavor os- 
cillation scheme where the predicted fluxes agree with the 
measurements in the GALLEX/SAGE, Homestake and (Su- 
per) Kamiokande experiments. The shaded regions and con- 
tours are the allowed regions at 95% CL. for different ax- 
ion losses (gio = shaded, 4.5 solid, 10 dashed and 15 
dash-dotted line). For the sake of clarity we omit the vac- 
uum solution for gio=15. 



ties of the seismic model. To constrain the axion-photon 
coupling through helioseismology we thus only take into 
account the differential effect of axion emission on solar 
models (Fig. ||). For g w = 10, the sound speed deviates 
from the reference model by 0.8% at the local maximum 
of the deviation at r ~ 0.1 Rq. This is probably too much 
to be absorbed by uncertainties in the input physics of 
the standard solar model. 

New measurements from MDI/SOHO providing more 
accurate low-1 frequencies should enable us to get a more 
precise determination of the sound speed in the central 
parts. But as the difference of the model with g 10 = 4.5 
to the reference model (gio — 0) is very small everywhere 
in the interior, these new frequencies should have almost 
no influence on our limits. 

Further, the <?io > 10 models produce less than 0.227 
for the surface helium abundance, while the smallest 
value allowed by helioseismology is 0.238. Put another 
way, we may no longer adjust the presolar helium abun- 
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FIG. 6. Expected rates for GALLEX/SAGE, Homestake 
and (Super)Kamiokande for different solutions and axion 
losses. Theoretical uncertainties are shown by la-error bars. 
For each g 10 < 10 we plot the SMA, LMA and VAC solu- 
tions from left to right; for gio = 15 the SMA and LMA cases 
have merged. The shaded bands represent the lcr-range of 
the measured rates. 



dance to any desired value to absorb the effect of axion 
losses on the solar model because the modified helium 
abundance is not compatible with seismological proper- 
ties of the convection layers. 

Finally, solar models with gio > 10 produce such large 
neutrino fluxes that the measured values cannot plausibly 
be explained by neutrino oscillations. While oscillations 
are not yet fully established as the undisputed expla- 
nation of the measured solar neutrino flux deficits, it is 
certainly not true that any neutrino flux can be absorbed 
by suitably adjusted mixing parameters. 

In summary, we believe that our present knowledge of 
the Sun excludes an axion-photon coupling in excess of 
gio = 10, corresponding to L a = 0.20 Lq, even if we can- 
not assign a statistically meaningful confidence level to 
this bound. This improves the previous solar limit 
of gio ^ 25 by about a factor of 3. While other exotic 
energy-loss channels would have a different energy and 
temperature dependence, it is probably generic that he- 
lioseismology excludes any such channel much in excess 
of the 0.2 L Q level. 

A recent helioscope experiment to search for solar ax- 



ions has reported a limit of gio < 6, valid for m a < 

0. 03 eV 0. Therefore, in this mass range the helioscope 
limit improves our bound and thus is self-consistent, 

1. e. its validity does not depend on an axion luminosity 
in excess of what is allowed by the properties of the Sun. 

Our new helioseismological bound is much weaker than 
the well-known limit 510 ^ 0.6 which has been derived 
from globular-cluster stars jl],|| . We still think it is useful 
to have independent information from the Sun, especially 
as the experimental accuracy of solar observations is or- 
ders of magnitude better than for every other star. Most 
importantly, the Sun serves as a source for experimental 
axion searches. Note that a proposed helioscope exper- 
iment using a decommissioned LHC test magnet could 
conceivably reach the globular-cluster limit Q . 
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APPENDIX A: ENERGY-LOSS RATE 

The axionic energy-loss rate of a nondegenerate plasma 
by the Primakoff effect was calculated in Ref. |32). We 
follow the representation given in Ref. [|l|, 



-in 



= 0.892 x 10- 3 erg g" 1 s" 1 g\ T% ft 1 F(k 2 ), (Al) 

where Tj = T/10 7 K and pi = p/10 2 g cm -3 . Screening 
effects are described by the dimensionless function 



F ( K ) = V^ dx ~ 
2?H J e 



x 



1 

2 1 .,1 



(x z 4- k ) In I 1 4 



(A2) 



which is shown in Fig. [?]. Here, 
n B 



16.56p 2 T 7 - 3 (V e +^^ 



(A3) 



where a = 1/137 is the fine-structure constant, ub the 
baryon density, Y e the number of electrons per baryon, 
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Yj the number of nuclear species j per baryon and Zj 
its charge number. In a standard solar model k 2 = 12 
throughout the Sun with a variation of less than 15%, 
corresponding to F(12) — 1.842. A good fit in this region 
is given by 



F{k 2 ) = 1.842 (k 2 /12)°- 31 , 
which is shown as a dotted line in Fig. 0. 



(A4) 
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FIG. 7. Function F(k 2 ), solid line exact result according 



to Eq. HAS]), dotted line our approximation of Eq. (A4) 
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